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Harper’s equation (aka the “almost Mathieu” equation) famously describes the quantum dynamics of an 
electron on a one dimensional lattice in the presence of an incommensurate potential with magnitude V and 
wave number Q. It has been proven that all states are delocalized if V is less than a critical value Vc = 2t and 
localized if k > Y. Here, we show that this result (while correct) is highly misleading, at least in the small Q 
limit. In particular, for V < Vc there is an abrupt crossover akin to a mobility edge at an energy E^', states with 
energy |£| < Ec are robustly delocalized, but those in the tails of the density of states, with \E\ > E^, form a 
set of narrow bands with exponentially small bandwidths ~ t exp[-(2;rQ'/2)] (where a is an energy dependent 
number of order 1) separated by band-gaps ~ tQ. Thus, the states with \E\ > E^ are “almost localized” in that 
they have an exponentially large effective mass and are easily localized by small perturbations. We establish 
this both using exact numerical solution of the problem, and by exploiting the well known fact that the same 
eigenvalue problem arises in the Hofstadter problem of an electron moving on a 2D lattice in the presence of a 
magnetic field, B = QI2n. From the 2D perspective, the almost localized states are simply the Landau levels 
associated with semiclassical precession around closed contours of constant quasiparticle energy; that they are 
not truly localized reflects an extremely subtle form of magnetic breakdown. 


I. INTRODUCTION 

Quasiperiodic ordered states of matter possess two or more 
periodic structures whose periods are incommensurate with 
each other. Examples of such materials include quasicrys¬ 
tals and crystals with incommensurate charge or spin den¬ 
sity waves as well as simple crystals in the presence of 
a generic uniform magnetic field. The absence of even 
discrete translational symmetry fundamentally distinguishes 
quasiperiodic from crystalline (periodic) structures. The ab¬ 
sence of quenched randomness makes them conceptually dis¬ 
tinct from disordered systems, although because in both cases 
Bloch’s theorem does not apply, there is no requirement that 
the elementary excitations are delocalized. 

As one of the simplest examples in the family, Aubry 
and Andre studied in 1980 the following one-dimensional 
finite-difference Schrodinger equation with a quasiperiodic 
potential^ 

t {fn+l + fn-l) + V COS (Qn + 0)f„ = Ef„ (1) 

where « 6 Z labels discrete sites, f is the nearest neighbor 
hopping amplitude, V and Q are the amplitude and wave vec¬ 
tor of the incommensurate potential (Ql'ln e (0,1) is irra¬ 
tional), and -6 /Q defines an “origin” of the incommensurate 
potential. Without loss of generality, we use the convention 
in which V > 0, f > 0 and open boundary conditions for a 
finite size system. This eigenvalue equation is generally re¬ 
ferred to in the physics literature as Harper’s equation or the 
Aubry-Andre equation and in the mathematics literature as the 
“almost Mathieu equation” (an analogy with the continuous 
Mathieu equation). We will use the name Harper’s equation 
from here on. 

Clearly, this is the eigenvalue equation corresponding to a 
one-dimensional tight-binding Hamiltonian in the presence of 
a quasiperiodic potential, which in second-quantized repre¬ 
sentation is 

^ ^ + S) clcn] ( 2 ) 


where E is an energy eigenvalue and - Yjfnc]i is the ere- 

n 

ation operator for the corresponding eigenstate. The same 
eigenvalue equation arises in the “Hofstadter problem” of 
an electron moving in a two-dimensional square lattice in 
the presence of a magnetic field corresponding to Ql'ln flux 
quanta per plaquette. 

There is a well-known self-duality in Harper’s equation^— 
which interchanges the roles of the kinetic and potential en¬ 
ergy terms, Vjl <-> t. Specifically, in terms of Fourier trans¬ 
formed variables. 
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where L is the system size. It is easy to see that satisfies 
the dual Harper’s equation 

U 

I^igm+i + gm-i) + It COS (Qm + ip)g,„ ^ Eg^. (4) 

This equation is self-dual when V - It. 

The spectrum and transport properties of this system have 
been extensively studied for over thirty years, in a wide set of 
contexts in mathematics and statistical mechanics. In partic¬ 
ular, the solution of the “Ten Martini Problem”ii^ proves that 
the spectrum is a Cantor set for V 4^ 0 and arbitrary “incom¬ 
mensurate” Q, i.e. so long as Qjln is irrational. There has 
also been progress towards experimental realization in cold 
atom systems^. It is widely believediriii^ that for any Q that 
is incommensurate, there is a single “metal-insulator” transi¬ 
tion at the self-dual point, V = It, such that the spectrum is 
absolutely continuous and all eigenstates are de-localized for 
V < It, while the spectrum is pure-point and all eigenstates 
are localized for V > It. The same line of reasoning leads to 
the conclusion that the spectrum is singularly continuous and 
the eigenstates are ‘critical’ for V - It. 
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In this paper, we revisit Harper’s equation from various per¬ 
spectives. We show that a conventional semiclassical analysis 
of the associated Hofstadter problem suggests that for a weak 
potential, V < 2t, and small QI2n <s 1, there exist two mobil¬ 
ity edges at £’c,± = + |2f - V\ that separate the localized states 
with E < Ec- and E > Ec,+, from the delocalized states with 
Ec- < E < Ec,+ ', i.e., the states near the band edges are ef¬ 
fectively localized even for a weak incommensurate potential. 
This conclusion is apparently verified by “exact” numerical 
studies, even when eigenenergies are computed with an accu¬ 
racy exceeding 1 part in lO'^. The resolution of this appar¬ 
ent contradiction comes from a calculable breakdown of the 
semiclassical analysis, and a careful analysis of the numerics; 
while the states outside the apparent mobility edges are, in 
fact, delocalized, they form exponentially narrow bands - es¬ 
sentially Landau levels - with correspondingly large effective 
masses, m* ~ r' exp[(2Q';7r/0], where a is a dimensionless 
function of \E — Ec\ and V72f. While these states are techni¬ 
cally delocalized, for even moderately small QI2ji they are so 
weakly dispersing that for all physical purposes they behave as 
if localized; we refer to these as “almost localized” states, un¬ 
derstanding the physics of which is one of the main results of 
the paper. Moreover, a similar analysis in the case V > 2t, still 
with QI2ji ^ 1, reveals that there is a hidden “almost transi¬ 
tion” characterized by the emergence of a real-space Fermi 
surface. For large wave vector QI2ji ~ (9(1) and relatively 
small potential V, we find that perturbation theory gives a sat¬ 
isfactory characterization and results in a series of gaps in the 
spectrum, whose sizes and locations are determined by the 
values of Vjt and Q\ all states are robustly de-localized. 

The rest of the paper is organized as follows. In Sec. m we 
study the physical properties of Harper’s equation in the limit 
where Q <& 2ti, analytically using a semiclassical approxima¬ 
tion and by numerical methods, which suggests the existence 
of two mobility edges. In particular, because the system is 
one dimensional, we are able to employ an extremely efficient 
recursive method to obtain numerical solutions for extremely 
large system sizes, even with lengths in excess of 10^ sites, so 
we have been able to test the validity of all our analytic argu¬ 
ments with great precision. In Sec. Hill we consider high-order 
perturbations and show that it reconciles the conflict between 
Sec. inland previous phase diagram. Especially, we character¬ 
ize the essential properties of the states near the band edges 
and why the conclusions of Sec. HJis physically relevant. In 
Sec. HVlwe discuss the previously overlooked “transition” as 
a function of energy that occurs when Q ^ 2n and Y > 2f, 
even though all the states are robustly localized. In Sec. lYl 
we briefly present our theoretical understanding and numeri¬ 
cal results for the situation in which the wave vector Q ~ (9(1). 
We conclude and discuss relations with previous work in Sec. 

lYll 


II. LONG PERIOD INCOMMENSURATE ORDER, Q « In 

In this section we present semiclassical and numerical 
studies concerning the spectral and transport properties of 
Harper’s equation when the period of the incommensurate po¬ 



Open along kx 

Closed 

Open along ky 

V>2t 

\E\<V- 2l 

\E\>V- 2t 

N.A. 

V <2t 

N.A. 

\E\ > 2t-V 

lEi < 2f - V 

V = 2t 

£ = 0 

EtO 

£ = 0 


TABLE I: Topologies of the constant energy contours. See Fig. [T]for 
illustrations. Note that for V = 2f and E - 0 the contour is a perfect 
square and connected along both the and ky directions at isolated 
points {n, 0) and (0, n). 


tential is large in units of the lattice constant. 


A. Semiclassical theory of a 2D crystal in a magnetic field 


It is well known that Harper’s equation also describes the 
Hofstadter problem - a model of a two-dimensional (2D) crys¬ 
tal in the presence of an incommensurate magnetic flux den¬ 
sity. To establish such an equivalence, consider a 2D tight- 
binding model on a square lattice with hopping matrix element 
t along the x direction and V/2 along the y direction and with 
an effective magnetic field B — Bz where B - QI2n so that 
there are QI2n magnetic flux quanta per plaquette. In a gauge 
chosen to preserve translation symmetry along the y direction, 
so that the Bloch wave-number ky is a conserved quantity, the 
Hamiltonian is 

Hid = X (^l+i + H-C. j H- V cos (fin H- ky'j clj^ Cn,k, 

n,ky 

(5) 

where n is the lattice distance along the x direction. Thus, in 
the subspace of fixed ky, with the identification 0 <-> ky and 
leaving the ky index implicit ^ c„, the Hofstadter prob¬ 
lem problem is seen to be identical to Harper’s equation in 
Eq. 12 Moreover, in the thermodynamic limit, when QI2n 
is irrational, spatially averaged physical quantities such as the 
density of states (DOS) and the localization length are inde¬ 
pendent of ky for the same reason that they are independent of 
the origin of coordinates (6) of an incommensurate potential. 
(See Ref. 0 and Appendix lAl for more detailed discussions.) 
Thus, the summation over ky in Eq. |5lonly produces a degen¬ 
eracy factor of Ly. 

Alternatively, we can choose the gauge that preserves trans¬ 
lation symmetry along the x direction; 


Hid’ 


Xi [y 2fcos(ki - Qm)cl i^c„t^ic, 

m,kx 


( 6 ) 


where m and kx label the lattice site and momentum along 
the y and x directions, respectively. In this gauge, E[ 2 D' is 
equivalent to the dual form of Harper’s equation in Eq. 12 The 
duality is thus identified with two different gauge choices for 
the Hofstadter problem. 

We can study the localization properties of the one¬ 
dimensional crystal with an incommensurate potential by con¬ 
sidering whether the eigenstates of the corresponding two- 
dimensional crystal in a magnetic field are localized along 
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the X direction, which is the spatial direction in the origi¬ 
nal one-dimensional problemiS. When the “field strength” 
B = Ql'ln <s; 1, the dynamics of Bloch electrons is given 
semiclassically by the Lorentz force law: 


elk . -> e 

fi— - -ev(k) X B - - 
at n 


Bx 


dE{k) 
dk , 


(7) 


where v(k) - is the electron group velocity. The re- 

fidk 

suiting electron orbits in two-dimensional momentum space 
are confined to constant energy contours determined by the 
zero-field energy dispersion E — eu - Itcosk^ + Vcosky. 
Since the velocity of an electron is perpendicular to the con¬ 
stant energy contour, the semiclassical motion is localized in 
the X direction unless the contour is open along the ky direc¬ 
tion, which occurs when V < 2t - lEI. (See Table Ufor the 
parameter regimes for different topologies of the Fermi sur¬ 
face and Fig. [T]for illustrations.) This seemingly implies that 
for potential amplitude V < 2t there are two mobility edges 
at £’c,± = ± (2f - y); the states are de-localized between these 
energies, while all states beyond the mobility edges are lo¬ 
calized. (For V > 2t the same analysis leads to the conclu¬ 
sion that all states are localized.) If we start from a pristine 
one-dimensional crystal and slowly turn on V, the states near 
the band edges localize first, and as V increases the energy 
window of de-localization between the two putative mobility 
edges becomes smaller and eventually closes at V = 2t. This 
is in sharp contrast to the expectation that all states are de¬ 
localized for V < 2t. 

The wave vector Q in Harper’s equation plays the role of 
a magnetic field in the two-dimensional crystal scenario. For 
B - QI2n <K 1, the semiclassical results are expected to be 
asymptotically exact, and the transition should thus be almost 
independent of Q. (Exceptions can arise when V «c 2f or 
y » 2f where the strong anisotropy and curvature of the 2D 
constant energy contours invalidates the semiclassical approx¬ 
imation, as will be discussed further in Sec. |III]) For larger 
Q, however, the magnetic field can induce tunneling between 
separate semiclassical orbits, a phenomenon known as mag¬ 
netic breakdownii. In particular, for Q ~ (9(1) there is no 
reason to trust the semiclassical theory at all, so for the pur¬ 
poses of the present section, we limit ourselves to small Q and 
will only return to consider 2 ~ (9(1) in Sec. 

To test the validity of our semiclassical arguments, we will 
now discuss numerical studies of the localization length and 
the density of states (DOS). 


B. Numerical results 





Since the one-dimensional Hamiltonian can be written as a 
tri-diagonal matrix, numerical calculations are efficient with 
recursive methods even for large system sizes^S. In particu¬ 
lar, we calculate the two-point Green’s function G («,«') = 
{E + i6-H)^\,. If the states are localized at energy E, 
then the Green’s function should be exponentially decay¬ 
ing with the distance between the two points G(n,n') oc 


FIG. 1: Constant energy contours (a) for V < It, (b) at the self¬ 
dual point, y = 2f, and (c) for V > 2t. The black curves are at 
E < -|2f- VI, the red at £ = 0, and blue at £ > |2?- V\. The 
black and blue curves are closed contours while the red one is open 
in different directions for different V and is at the Lifshitz transition 
when V = 2t. Note that ( = 1.0, V = 2.5 and t - 0.8, y = 2 are 
equivalent up to an energy rescaling: £ —> 0.8£. 
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E 

FIG. 2: The inverse localization length l/T as a function of energy 
for three different parameter regimes at Q = 1/31. The system size 
is L = 7.2 X 10^. The blue and black curves are related by the duality 
transformation in Eq. [3]connecting the Hamiltonians in Eqs. [T]and 
|4] The localization length is always finite for the blue curve, for 
which V > 2t and corresponds to the constant energy contours in Fig. 
|l(c)| However, there is a portion of the spectrum that is apparently 
delocalized for the black curve, for which V < 2t corresponding to 
the constant energy contours in Fig. |l(a)| The red curve at the self¬ 
dual point V = 2l corresponds to the contours in Fig. |l(b)| 


exp(- \n - n'\ jA) where A is the localization length; other¬ 
wise the system is delocalized and T ^ oo (l/T ^ 0). In 
all cases we have carried out calculations on systems of sizes 
L large enough that the results are independent of L - which 
often means up to sizes L ~ 10® sites. 

We study the localization properties for a representative 
small Q — 1/31 = 27r(l/627r) ~ 27r(5.1 x 10“^) and for var¬ 
ious values of f, V, and E that lead to different topologies of 
the associated two-dimensional constant energy contours as 
shown in Fig. [1] In Fig. I^we display the results of numer¬ 
ical studies of the inverse localization length as a function of 
E for various values of Vjt. The results are consistent with 
the inferences made in the previous section on the basis of the 
two-dimensional semiclassical theory; if the energy contour 
is closed (black and blue curves in Figs. |l(a)||l(c)] l or open 
along the direction (red in Fig. |l(c)[ ), the corresponding 
states in the one-dimensional (ID) problem are localized. If 
the energy contour is open along the ky direction (red in Fig. 
|l(a)| i, the ID states are delocalized. In particular, for V < 2t 
but IFI > 2t - V the results indicate that the states are local¬ 
ized - consistent with expectations from semiclassical theory 
but in conflict with previous claimsk^i^ii. 

We have also computed the density of states (DOS) at en¬ 
ergy E, defined as p{E) = using the same 

n 

methodsiS. Our algorithm is most efficient and accurate when 
the system is localized and the Green’s function is exponen¬ 
tially suppressed. However, even if the parameters to be con¬ 
sidered are in the de-localized regime, the spectral properties 
are the same as those of the dual Harper’s equation given in 
Eq. m which is necessarily localizedi^. For this reason, where 



E 



E E 

FIG. 3: The DOS p{E) of Harper’s equation with parameters r = 1, 
F = 1.6 and Q = 1/31. The system size is L = 7.2 x 10^ and 
6 = 0.0001 brings in a small but finite energy resolution. Drastically 
dilferent behaviors are observed above and below E = ±0.4. Note 
that the spectrum is particle-hole symmetric therefore only the £ < 0 
half is shown. In practice, it is obtained through the computation for 
the dual Harper’s equation with t = 0.8, V = 2 and Q = 1/31. The 
lower panels are regimes between and just beyond the apparent mo¬ 
bility edges on an expanded scale, another expansion near the band 
edge is shown in Fig. 


the states are delocalized, we always compute the DOS from 
the localized dual. 

An example of the sort of DOS so obtained is shown in 
Fig. |3]for f = l,y = 1.6 and Q - 1/31. Drastically dif¬ 
ferent behaviors are observed above and below E - ±0.4. 
For IFI > 0.4, there are discrete sharp peaks separated by re¬ 
solvable gaps while for 0.4 > IE], the DOS appears to be a 
smoothly varying function of E. The correspondence with the 
semiclassical treatment of the 2D version is that the discrete 
peaks correspond to Landau levels where the constant energy 
contours (at zero field) are closed. The continuous spectrum 
arises where the energy contours are open, and the apparent 
mobility edges correspond to the Lifshitz transitions at which 
the topology of these contours changes (see Table|I]l. 

More generally. Fig. |4] is the phase diagram of Harper’s 
equation with Q - 1/31 extracted from localization length 
calculations. The phase boundaries for V <2t and V > 2f are 
colored in black and red, respectively, and are fully consistent 
with Ec - ± (2t — V) (blue lines) derived from the semiclas¬ 
sical treatment of the 2D model, as given in Table |I] On the 
V <2t side of the phase diagram, an apparent metal insulator 
transition produces qualitative changes in both the localiza¬ 
tion length A and the DOS p. On the V > 2t side, however. 
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FIG. 4: The phase diagram for Harper’s equation with f = 1 and 
2 = 1/31 obtained from numerical study of localization length on a 
system of size L = 7.2 x 10^. For each V < 2t, Ec (black) separates 
the regimes of metal (localization length A oo) and insulator (i.e. 
with available numerical accuracy, A is hnite). For V > 2, the criti¬ 
cal values of (red) are obtained by a duality transformation of the 
V <2t cases. Change of DOS behaviors is observed on both V < 2t 
and V > 2t phase boundaries, though on both sides of the latter phase 
boundary the system is localized and insulating. The blue solid lines 
are = ± (2r - V) in accord with Table |I] The orange dotted lines 
are the band edges at £ = ± (2f -l- V). One important caveat is that 
there is an implicit energy resolution on the level of machine preci¬ 
sion in our numerical calculations. For infinite precision, the phase 
diagram will essentially reduce to a single metal-insulator transition 
at V = 2t, see Sec. |nI]for more details. 

while the analogous change in the structure of the DOS seem¬ 
ingly signifies the existence of a transition (the dual of the 
metal-insulator transition) the localization length is finite on 
both sides of the transition. We will discuss the nature of this 
“insulator-insulator” transition in Sec. 

C. Localization near the band edges from the 2D perspective 

It is easy to prove that the spectrum of Harper’s equation is 
confined to the range of energies lisl <2t+V. One important 
qualitative point that is less obvious is that the states near the 

lower band edge E - 2t - V (and the upper band edge, E ~ 

2t+V) always appear localized, whether or not V is larger than 
2t. This is most easily understood from the 2D perspective, 
H 2 d in Eq. 0 In the absence of an effective magnetic field 
{QI2n - 0), the dispersion near the band-bottom is accurately 
treated in the effective mass approximation, which means that 
for small but non-zero Q, the low energy spectrum is well 
approximated as Landau levels. This accounts both for the 
discrete peaks in the spectrum and their spatial localization. 

To be specific, we expand the zero-field energy dispersion 
of H 2 d to quadratic order of k around the band minimum at 
(7r,7r); 

, V 2 

ek - 2t cos kx + V cos ky ~ -2t - V + tq^ + —q^ + ... (8) 


— Numerics 
I Quadratic order 



E 


FIG. 5: The DOS p(£) with the same parameters t = I, V = 1.6, 
6 = 0.0001, and Q = 1/31 as in Fig. [3 near the lower band edge 
but on an expanded scale.. The red and blue marks on the top label 
the locations of the discrete Landau levels from the expansions of the 
electron dispersion to quadratic order (free-electron approximation) 
and quartic order in momentum, respectively. 


where (qx, qy) - (kx, ky) - (n, n). The effective magnetic field 
is QI2n flux quanta per plaquette, and the corresponding ef¬ 
fective cyclotron frequency is uc - Q y/2tV. Thus, the energy 
levels near the lower band edge are quantized in Landau levels 
with e„ ^ -2t -V + (n+ \l2)uc where h = 0, 1, 2, ■ ■ ■. We 
show in Lig. |5]an enlargement of the DOS near the lower band 
edge obtained from numerical solution of //lo with f = 1, 
V - 1.6, and g = 1 /31 as in Lig. [3j the corresponding values 
of e„ are indicated by the red marks at the top of the figure. 

The consistency is remarkable for the few lowest energy 
levels, where the effective mass approximation is highly ac¬ 
curate; at somewhat higher energy, however, while the en¬ 
ergy levels remain sharp and gapped, the level spacing grad¬ 
ually decreases. However, the accuracy of the theoretical es¬ 
timates can be improved by taking into account the higher or¬ 
der terms in the expansion of the energy dispersion, e^, e.g. 
V cosky ~ - V + jqy - ^qy + 0(qy). The Landau level prob¬ 
lem including the quartic terms is discussed in Appendix |B] 
and the results (obtained with minimal numerical worbiS) are 
shown as the blue marks at the top of Lig. |5] This accounts 
quite accurately for the positions of the peaks in the DOS 
over the entire range of energies shown in the figure. The 
spatial extent of the wave-functions is determined by the ef¬ 
fective magnetic length, { oc which also determines 

the relevant range of momenta Ak ~ Ijt. Thus, the valid¬ 
ity of the expansion in powers of k, and hence the validity 
of the Landau level analysis, relies both on the smallness of 
Q and on E being close to the band edge, wc ^ 4f and 
\E + 2t + V\ <s At. However, at this level of approximation, no 
aspect of the analysis is sensitive to whether or not V is larger 
or smaller than 2t (as long as the Landau level spacing uc 
is much smaller than the original band width, which requires 
(QI2ti)^ «c y/2f «c (2nlQ)^). We will focus more carefully 
on the properties of these band-edge states in Sec. m 
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Closer to the band center, the fact that there is a periodic lat¬ 
tice qualitatively affects the Hofstadter problem. As required 
by duality, the resulting change from closed to open topology 
of the constant energy contours that occurs at £ = ±Ec pro¬ 
duces similar changes in the character of the DOS for V < 2t 
or y > 2t. For V < 2t, there appears to be a continuous DOS 
for l^l < Ec and the states are delocalized in the x direction, 
so in terms of the properties of the original ID problem, Ec 
acts as a mobility edge in the traditional sense. However, for 
V > 2t, from the perspective of the original ID problem, the 
states are robustly localized both above and below Ec', the na¬ 
ture of the crossover that occurs at £ = ±Ec in this case will 
be discussed in Sec. m 


III. RESOLVING THE PARADOX - DELOCALIZATION 
OE THE STATES NEAR THE BAND EDGES 

As we have shown in Sec. [Ill the semiclassical and nu¬ 
merical results disagree qualitatively with the previous litera¬ 
ture concerning both the localization and spectral properties of 
Harper’s equation for small Q <& 2n. In this section, we first 
briefly summarize the previous arguments that there is a sin¬ 
gle energy independent metal-insulator transition at V = 2t, 
and then examine more carefully the properties of the states 
near the band edges where the contradiction arises. In partic¬ 
ular, we focus on two perturbations that have been neglected 
in the semiclassical approximation and show how they affect 
the localization problem. We will see, for Q^2n and Y < 2t, 
the states near the band edges behave ‘almost localized’; even 
though their wavefunctions are extended, their bandwidth are 
exponentially small in I IQ, therefore they are easily localized 
in the presence of small perturbation or finite energy resolu¬ 
tion. 


A. Previous arguments on the localization transition 

First of all, we briefly summarize the previous argumentsii^. 
In the notation of Eq. [3] we define 

fix) = exp iiniQx + iOm) (9) 

as an extension of fn, n 6 Z to .r 6 R. Since m takes integer 
values, for a g„, which remains normalizable in the L ^ oo 
limit (that is, g is localized), f{x) is a Bloch function with 
a fundamental period of 2nlQ. Then the wavefunction /„, 
though no longer periodic, is still extended due to its periodic 
envelope/(x). Namely, if the solution to one specific Harper’s 
equation is localized, its dual is necessarily extended. 

We can further use the fact that the inverse localization 
length is related to the DOS through the Thouless formula: 

l/d(£) = J dE'lnj(E-E')/tjp(E') (10) 

which can be equally applied to the dual system in Eq. |4]with 


the same DOS p: 

1/T'(£) = J ^/£'ln|2(£-£')/y|p(£') (11) 

then !//!(£) = 1/A'(E) + ln(y/2f). This relation holds well as 
can be seen in Eig. [J] the Harper’s equations corresponding 
to the blue and black curves are dual, and indeed there is a 
constant vertical shift of ln(y/2f) between them. 

Let us first consider the case where states in the original 
system are localized: its dual is extended so !//!'(£) = 0, 
therefore 1/A(E) - ln(y/2f) > 0 is constrained to y > 2f. 
Similarly, if the original system is extended !//!(£) = 0, its 
dual !//!'(£) = ln(2f/y) > 0 requires V < 2t. Consequently, 
all states are localized at all energy levels £ for V > 2t, and 
all states are extended at all energy levels £ for V < 2t. 

There is an obvious contradiction between this argument 
and the two-dimensional semiclassical problem, where the du¬ 
ality corresponds to switching the roles of x and y. For open 
constant energy contours (red and blue curves in Figs. |l(a)| 
and |l(c)[ l, the semiclassical theory does suggest that the sys¬ 
tem is either localized along x and de-localized along y, or lo¬ 
calized along y and de-localized along x. However, there is a 
third possibility that the constant energy contour can be closed 
(black curves in Figs. |l(a)| and |l(c)| l, where the states are 
physically localized along both the x and y directions. More 
specifically, while the previous claim that all states are local¬ 
ized for V > 2t and the states between the two apparent mo¬ 
bility edges, i.e. those with Ec = +(2f - y) are de-localized 
for y < 2f is consistent with the results in Sec. im both the 
semiclassical theory and the numerics seemingly imply that 
the states beyond the mobility edges at +£ e (2r - y 2f -i- y) 
are localized even for V < 2t. 


B. Localization near the band edges from a real-space 
perspective 

In Sec. Ill Cl and more generally in Sec. Ill Al we found 
that the states near the lower band edge form effective Landau 
levels; at this level of approximation, the states in a Landau- 
level are exactly degenerate, and thus it is possible to construct 
eigenstates that are either localized or delocalized. In this cir¬ 
cumstance, even parametrically small terms that are omitted 
can, in principle, play a qualitative role in resolving this de¬ 
generacy. In this subsection we consider the states near the 
lower band edge from a ID real-space perspective, treating 
2 «c 2;7r as a small parameter. In particular, we show that 
exponentially small terms, ~ exp[-(2a;7r/2)] lift this degen¬ 
eracy in such a way that the band-edge states are delocalized 
for V < 2t and localized for V > 2t. Note that terms that pro¬ 
duce these behaviors are so small for even moderate values 
of 27t/Q that they are entirely unobservable in the numerical 
studies of Sec. IIIBl despite their being carried out to less than 
1 part in 10^^ uncertainty! Conversely, since it is difficult to 
imagine circumstances in which the energy uncertainty (due 
to thermal broadening if nothing else) is anywhere near this 
small, this observation carries with it the implication that/or 
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all plausible physical purposes, these states are ejfectively lo¬ 
calized. It is in this sense that they are “almost localized.” 

Consider the states within a quantized energy level e„ and 
the corresponding Hilbert subspace consisting of a single lo¬ 
calized state per potential well consisting of 2nlQ sites. 
There are two perturbations that are exponentially small in 
1/2 and implicitly neglected in the semiclassical theory in 
Sec. M- the tunneling between the potential wells t^fi and 
the local energy variation between the potential wells Veff, 
which turn out to be the key to resolving the conflicting con¬ 
clusions. The way this occurs can be most easily seen from 
a simple variational analysis of the band-bottom states in the 
limit 2 « 1. 

Local energy difference between wells: To begin with, we 
construct a variational state centered near a single minimum 
of V of the gaussian form 

/„ = yve™ exp[-(u - (i2) 

where N is the normalization constant, while the center of lo¬ 
calization, a, and the localization length, f, are treated as vari¬ 
ational parameters, although we will assume (and then con¬ 
firm) that 2nlQ'3> ( :s> The variational energy in this state 
is easily seen to be (using the Poisson summation formula) 

{H) = Vcos{Qa)e-^"'^"^'^ - 2te-^l'^' + 6E (13) 

'> t c2 

where 6E is all terms of order with integer in > 1. 

Clearly, the leading order terms (ignoring 6E) simply repro¬ 
duce the results of the harmonic approximation in Eq. [8] The 
result is a = ;r(l +2m)IQ,( - (2f/y)*^"^-yT/S giving a varia¬ 
tional energy eo ~ -2-t - V -i- 2 ^2tV /2 that is independent of 
m, i.e. the minimum of the potential in which it is localized. 

This result is not exact. The leading contribution to its cor¬ 
rections can be readily seen by looking at the leading m de¬ 
pendent term in 6E\ 6E - -2ye^’''^ cos(27rfl) H-... where ... 
refers to other terms of the same order which are independent 
of a, as well as higher order terms in powers of . Thus, 
we see that there is an exponentially small difference in energy 

V,fi ~ V exp[-7T^f] (14) 

between a state localized near a — njQ and the neighboring 
well at a = 37r/2- 

Tunneling between the potential wells: Using the same vari¬ 
ational wavefunctions, we can readily estimate the tunnel¬ 
ing matrix element between two neighboring states by sim¬ 
ply evaluating the matrix element of the tunneling term in the 
Hamiltonian between localized states centered at a - njQ and 
a = 37r/2- The result is 

fefF~fexp[-7rV2"^"]- (15) 

this is qualitatively consistent with the more careful analysis 
included in Appendix ICl 

From this simple variational treatment, we conclude that 
Ueff » left {he. the states are localized) so long as 2^^ < 1> 
and conversely, that Veff «c feft if 2^^^ > 1- Miraculously 
(and probably accidentally), since this variational approach 


yields - 2tfV, it reproduces the exact criterion for lo¬ 
calization. (Note that in terms of parametric dependence on 
2, ~ 2 which implies that both and are small 

compared to 1 when 2 ^ 1> promised.) 


C. Localization near the band edge from the perspective of 
magnetic breakdown 

To locate the localization transition more accurately, we 
consider related corrections to the 2D semiclassical theory 
from Sec. Ill Al So long as the effective magnetic field, QI2n, 
is non-zero, the semiclassical theory is not exact, and in partic¬ 
ular there is always a non-zero amplitude for magnetic break¬ 
down across the Brillouin zone boundaries. While the ampli¬ 
tude for such processes is exponentially small, they are exactly 
the terms that determine whether the states in a given Landau 
level are localized or not. 

Magnetic breakdown along the y and x directions produce 
essentially different outcomes: the former results in an open 
constant energy contour along the y direction and hence de- 
localization along the x direction, and thus is associated with 
feff; the latter results in an open constant energy contour along 
the X direction and hence is associated with local energy dif¬ 
ferences in potential wells, Vefj. Manifestly, when V > 2t, the 
constant energy contour near the band edge is anisotropic and 
elongated along the x direction, and conversely for V < 2t. 
(See Fig. [T]) Therefore magnetic breakdown occurs predom¬ 
inantly along the x, so that the states are localized suggesting 
Veff > 2feff, or y direction so that the states are de-localized 
suggesting Veir < 2feff, depending on whether V 72 f is greater 
or less than 1. Indeed, the constant energy contours are sym¬ 
metric under n/2 rotation only at the self-dual point of V = 2f, 
which means that only here is VefF/2fefF = 1. A further impli¬ 
cation of this is that the spectrum forms a Cantor set, since for 
each Landau band of the original problem, i.e. for each value 
of the energy level e„, the effective Hamiltonian confined to 
the space of states spanned by this Landau level defines a new 
version of Harper’s equation, but with much smaller energy 
scale and longer length scale. These results are indeed consis¬ 
tent with previous conclusionsi^i^iiS. 


D. Numerical tests for “moderate” Q/2n 

To test this perspective, we numerically exact diagonalize 
Harper’s equation Hamiltonian for a representative value of 
2 = 1/5, which is still small compared to 27r, but not so 
small that exponentially small effects are beyond the realm 
of numerical studies (even by setting 6 - 0). We take 
L - 377 ~ 12(27r/2), the probability distribution of some 
lowest eigenstates for f = 1 and V - 1.6 (upper panel) or 

V — 2.5 (lower panel) are shown in Fig. |6] We can clearly 
see that the states near the lower band edge are localized for 

V - 2.5t > 2t but clearly extended for V - 1.6f < 2f - once 
we can resolve the exponentially small perturbations. 

We note that similar arguments can be useful for the identi¬ 
fication of the metal-insulator transition and mobility edges in 




8 



-Kl' 



FIG. 6: The probability amplitude of lowest eigenstates of Harper’s 
equation for Q = 1/5, f = 1, V = 1.6 (upper panel) and V = 2.5 
(lower panel). The results are obtained through exact diagonalization 
on system size L = 377 ~ 12 (InIQ) with almost 12 complete poten¬ 
tial wells. The dotted red curves are illustrations of the V cos (Qx) 
potential profile (not to scale). 


the generalized versions of Harper’s equations such as the ex¬ 
amples considered in Ref. where one can simply compare 
the quantum tunneling amplitudes across the x and y direc¬ 
tions in the equivalent two-dimensional Hofstadter’s problem. 
Further details and examples are included in AppendixlDl 


IV. HIDDEN CROSSOVER FOR V > 2t AND Q « 2n: A 
REAL-SPACE FERMI SURFACE 

In the previous sections, we have discussed the effective 
metal-insulator transition for Harper’s equation with V < 2t 
and Q <Si 2n, where drastically different behaviors of the spec¬ 
tral and transport properties are observed in the physical limit 
on the two sides of the apparent mobility edges ±Ec - 2t-V. 
In the current section, let us turn to the related ‘insulator- 
insulator’ transition at ±Ec - y-2f in the duality transformed 
systems for V > 2t and Q <& 2ti, see Fig. |4] All states are lo¬ 
calized for y > 2f, yet the behavior of the DOS still shows a 




FIG. 7: The local electron integrated DOS of Harper’s equation 
dE\f„(E)f with f = 1, y = 2.5 and Q = 1/31 for states 
within the energy range -0.8 < E < -0.7 (upper panel) and 
-0.4 < E < -0.3 (lower panel), respectively. The red dotted curve 
is an illustration of the incommensurate periodic potential (not to 
scale). The blue arrows in the lower panel mark additional dips ap¬ 
pearing at the potential minima. 

sharp change in character at E^, see Fig. [3] What is the nature 
of this transition? 

First of all, the Pauli exclusion principle forbids the elec¬ 
tron density per site exceeding one. As we have shown in Sec. 
Hn near the band edges the electron eigenstates of Harper’s 
equation localized in each potential well resembles that of a 
harmonic oscillator with their electron density mostly concen¬ 
trated around the potential minima. Heuristically, for a deep 
potential well, the electron density there will reach one elec¬ 
tron per site after filling a finite number of quantized levels. 
For Q<k2ti the potential around the minima is slow varying, 
therefore the nearby lattice sites have a band width approxi¬ 
mately equal to VF = 4f - the band width of an infinite ID 
chain with y = 0. From the lower band edge E - —2t — V, 
these sites will reach maximum electron density at £ -2t-V 
consistent with Ec \fV> 2t. 

For illustration, we calculate the local electron inte¬ 
grated DOS of Harper’s equation within an energy range 















































































































9 


|/„(£)|^ with f = 1, y = 2.5 and Q = 1/31, 
and the results are shown in Fig. [T] for the energy ranges 
-0.8 < E < -0.7 (upper panel) and -0.4 < E < -0.3 
(lower panel) corresponding to the two respective insulating 
phases in Fig. |4] Characteristic dips in the local electron den¬ 
sity at the potential minima are clearly present as marked by 
the blue arrows in the lower panel, a signal that the electron 
density there has reached maximum below E - -0.4. These 
dips broaden as energy E increases and more sites reach max¬ 
imum electron density. Consequently, for Fermi energy Ef 
between ±Ec, the real-space electron density after filling all 
states below Ef receives non-analytic singularities when it 
drops from maximally filled to partially filled - similar to the 
momentum-space electron density of a metal at the Fermi sur¬ 
face but in the real space instead. The real-space electron den¬ 
sity \fn{E)'^ for all states below the Fermi energy Ef 

is shown in Fig. [Hlforii/r = -0.7 (upper panel) andZif = -0.3 
(lower panel). 

In Sec. m we have shown that the effective metal-insulator 
transition for V < 2f is accompanied by a qualitative change 
of DOS profile as the influence of the lattice comes into play. 
Duality suggested that these should apply equally to the tran¬ 
sitions at y > 2t. We note that our arguments in Sec. |II] 
near the band edges hold irrespective of V, therefore the band- 
edge states have the same properties no matter V > 2f or 
y < 2t: the eigenstates are effectively and independently lo¬ 
calized within each potential well regardless of the lattices, 
therefore the quantized and almost degenerate energy levels. 
However, once the electron density reaches maximum at the 
potential bottoms, further Ailing is gradually forced to the 
neighboring sites, which occurs sequentially in different po¬ 
tential wells according to their respective lattice displacement, 
changing the spread of the energy eigenvalues and the behav¬ 
ior of DOS. 

Note that this transition is unique for V > 2t. For V < 2t, 
the shallow potential means that the tunneling between the 
wells becomes relevant and the gaps in the local integrated 
DOS between the wells close before the maximum electron 
Ailing is ever reached anywhere. The transition is then dom¬ 
inated by the metal-insulator transition. Quite interestingly, 
the phase transition for V > 2f is an analogy of the metal- 
insulator transition for V < 2t but with the roles of real space 
and momentum space interchanged. For V < 2t and energy 
IFI <2t —V between the apparent mobility edges, the system 
is metallic in the sense that the real-space two-point Green’s 
function is long-range and the momentum-space electron den¬ 
sity is singular. On the other hand, for V > 2t and energy 
IFI < y - 2t between the critical points (the insulator* phase 
in the phase diagram in Fig. |4]i, the system has singular real- 
space density and the electron correlation is long-range in 
Fourier space. See Appendix |E] for more analytic details. We 
would like to emphasize that there is no symmetry breaking in 
either transition. Rather, they are identified by the appearance 
and disappearance of singularities and changes of physical ex¬ 
ponents, yet both are slightly rounded off in the presence of a 
finite incommensurate Q. Both transitions can be viewed ac¬ 
cording to the picture used in Sec. Ill Al as a Lifshitz transition 
in two dimensions projected to one dimension, but along dif- 




FIG. 8: The real-space electron density of Harper’s equation 
^dE If„{E)f with r = 1, y = 2.5 and Q = 1/31 for all states 
below the Fermi energy Ef = -0.7 (upper panel) and Ef = -0.3 
(lower panel), respectively. The blue dashed line indicates the maxi¬ 
mum electron density allowed - one electron per site. The red dotted 
curve is an illustration of the incommensurate periodic potential (not 
to scale). 

ferent directions. The duality transformation interchanges the 
two directions, therefore not surprisingly, it also exchanges 
the coordinates in which the transition occurs. 


V. LARGE Q LIMIT 

Up until this point, our discussion has focused on the limit 
in which 2 « 27r is relatively small. For large QI2n ~ (9(1), 
the effective magnetic field is large in the corresponding two- 
dimensional Hofstadter problem, so the semiclassical theory 
ceases to be useful. Conversely, we will consider the case in 
which y/2f is small, so in the original ID formulation of the 
problem, we can compute the effects of the incommensurate 
potential perturbatively and compare with numerical results. 
We And (consistent with our own earlier flndings in Ref. [l^ 
that 1) the lowest orders in the perturbation theory describe 
Harper’s equation satisfactorily so long as y/2f is even mod- 
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FIG. 9: The inverse localization length (upper panel) and the DOS 
(lower panel) for Flarper’s equation with r = 1, V = 0.2,0.5,0.6 and 
a large Q = 2^- (VS - l)/2 over the energy range E e [-2,0]. Due 
to particle-hole symmetry, only the E < 0 half is shown. The peaks 
in the DOS correspond to the van Flove singularities. 


erately small, and 2 ) all states are de-localized and form a 
continuous spectrum separated by the gaps, consistent with 
previous expectations. 

In the perturbation theory, we begin with a nearest-neighbor 
tight-binding model of a one-dimensional chain with energy 
dispersion - 2t cos{k), and then consider the gaps induced 
by the lowest order scattering from the incommensurate po¬ 
tential. More specifically, the nth order scattering generates 
an off-diagonal mass term between momentum k and k + iiQ, 
resulting in a gap of size ~ V" fW"^^ centering on energy 
En — - Sk+nQ - ±2cos{nQI2). The DOS vanishes in the 

gap and the localization length is finite; between the gaps the 
spectrum is continuous, and the localization length is infinite. 

To check these conclusions, we numerically compute the 
inverse localization length and DOS of Harper’s equation with 
a large Q. One example is shown in Fig. |9]with hopping am¬ 
plitude f = 1, potential amplitudes V - 0.2, 0.5, and 0.6 re¬ 
spectively, and wave vector (favored in much of the literature) 
Q - 27r ■ ( V5 - l)/ 2 . A series of gaps are clearly observed. 


and the inverse localization length is zero everywhere except 
inside the gaps. We note that the centers of the gaps from the 
first three orders of perturbation £1 = +2 cos(2/2) = +0.725, 
£2 = ±2cos(0 = +1.475 and £3 = +2cos(3e/2) = +1.793 
and their corresponding sizes are remarkably consistent with 
our numerical results in Fig. |9] In addition. Van Hove singu¬ 
larities in the DOS are clearly seen in the DOS on the band 
edges. These conclusions are straightforward to generalize to 
other Q. 

One major distinction between the incommensurate and 
commensurate potentials is that in the former case there is 
an unending cascade of higher order gaps while in the latter 
case new gaps are not generated beyond the order of perturba¬ 
tion theory corresponding to the commensurabilit}^. So why 
can we expect a finite-order perturbation theory to work for 
Harper’s equation? Heuristically, the obtained gaps as well 
as the resulting inverse localization length are exponentially 
small (~ (V72f)") which eventually results in gaps that are un¬ 
observable given a finite energy resolution 6 and system size 
L. In general, for any given level of desired energy resolution 
6 and system size L there is no essential distinction between 
the incommensurate problem and a suitable high order com¬ 
mensurate approximant. 


VI. CONCLUSION AND DISCUSSIONS 

We have revisited important properties of Harper’s equa¬ 
tion. This problem exhibits many extraordinary features, valid 
as matters of principle, if studied with infinite precision, i.e. 
assuming infinite energy resolution, perfect periodicity of the 
incommensurate potential, and no higher harmonics to the 
potential. This complexity includes conclusions concerning 
the non-existence of a true mobility edge which is a conse¬ 
quence of an exact self-duality, and a self-similar structure of 
the spectrum. However, if there is a small parameter, from a 
physical perspective, in which there is always a finite limit 
with which physical quantities can be controlled, exponen¬ 
tially small effects are physically negligible. 

For Q « 2n, the previous conclusion that there is a sin¬ 
gle metal insulator transition at V = 2 f for arbitrary energy £ 
holds in the strict sense. However, in the more physical sense 
in which a finite energy resolution is present due to disorder, 
finite temperature, system size limit etc., we find two “effec¬ 
tive” mobility edges at ±Ec — 2t—V for V <2t and Q <sc 27t; 
the states near the band edges are for most physical purposes 
localized even for a weak incommensurate potential. More¬ 
over, we have discovered a crossover at ±Ec — V — 2t for 
V > 2t and Q 2tt that is dual to the metal insulator transi¬ 
tion and characterized by the emergence of a real-space Fermi 
surface. 

Through both perturbative methods and equivalent mod¬ 
eling, we have self-consistently re-established the properties 
of the various phases of Harper’s equation. We note that for 
small Q <sc 2n, the two-dimensional Hofstadter problem has a 
quantized Hall conductance of cr^y = ne^ fh if we fill n discrete 
energy levels near the band edges, so that n electrons will be 
transported along the x direction in each cycle of ky e [0,27r]. 
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On the other hand, in the one-dimensional representation of 
Harper’s equation, ky is the initial phase of the incommen¬ 
surate potential, which is adiabatically shifted by exactly one 
period in each cycle of ky. Since there are n electrons localized 
in each period of the potential for the limit we are considering, 
it is straightforward to see that n electrons are pumped from 
one end of the chain to the other. Such adiabatic pumping 
process in quasiperiodic systems has been previously studied 
in Ref. [l2lfo Similar argument also implies a connection of 
Harper’s equation with certain commensurate wave vectors to 
three-dimensional Weyl semi-metalii. 

The methods and conclusions in the main text can be gener¬ 
alized straightforwardly to models with further neighbor hop¬ 
ping or more complex potentials^S. One interesting example 
is the ‘correlated disorder’ of the form V cos (Qx^) with an 
incommensurate wave vector Q. For y - 1, the potential is re¬ 
duced to Harper’s equation. Based on a different technique^^. 


mobility edges have been identified for y < 1 and the states 
are fully localized for y > 1. 
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Appendix A: The accuracy of 6 independence in Harper’s 
equation 

In this appendix, we discuss the applicability of physical 
quantities’ 0 independence. Mathematically, it was proven 
that the spectrum of the almost Mathieu equation (Harper’s 
equation) is independent of the choices of 0 in the thermody¬ 
namic limit and when QI2n is indeed irrationa l'll*^ . Here we 
offer a heuristic physical argument. 

To start with, we consider an intensive physical quantity 
Ai{0) for a one-dimensional system with finite size ( » 1, 
where exists in general a sizeable 0 dependence. For a 
larger system size L — nl, since the interfaces are negligi¬ 
bly small in comparison with the bulk, we may safely sepa¬ 
rate the measurement into various smaller components each 
with length I but different phases 0,„ - 0 + (m - l)2/(mod27r): 
Ai^ni{0) - therefore the 0 dependence is averaged 

m 

out over many different 0,„’s. We may then repeat the pro¬ 
cedure for even larger system sizes, and in each step the 0 
dependence scales down. This suggests that the 0 dependence 
is inverse proportionally suppressed as the total system size 
L. Similar argument works for extensive physical quantities 
Al=„i(0) - which receives just an additional factor 

m 

of the total system size L. 

It is straightforward to see that this argument is valid for the 
DOS p (per site), since the electron density (per site) is exten¬ 
sive (intensive). On the other hand, we show in the following 
that the correct self-averaging quantity for the localization is 
the inverse localization length l/T. 

The localization length is defined in terms of the Green’s 
function G by G{l,l;0) oc exp(-///i/(0)). Therefore, for a 
larger system L - nl: 

G(1,L;0) oc j~jG(m/— / H- 1, m/; 0) 

m 

ncd , /; 0„i) oc exp(-^(//I/(0m)) (Al) 


OC 
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FIG. 10: The system size effect of the spectrum dependence on 
different choices of 9 of the inverse localization length l/T (upper 
panel) and the DOS p (lower panel) of Harper’s equation with t = 1, 
V = 2.5, 2 = 1/31 and E = -0.3. The red dotted line is the average 
value over the spectrum for the largest system size. 


Comparing this with G(l,L;0) oc exp(—LIAi(6)), we have on 
the order 0(L): 


- L/Al{0) - -YjIAiiOm) 

m 

l/dL( 0 ) - - V l/d/( 0 ,„) 


(A2) 


We illustrate in Fig. [T0]the system size effect of the spec¬ 
trum dependence of the inverse localization length l/T and 
the DOS p on different choices of 9. For the system sizes we 
study in the main text, the 9 dependence has been suppressed 
well enough. 


mation is fairly consistent with numerical results on the DOS 
near the band edges. In this appendix, we show that the con¬ 
sistency can be made better over a wider range of energy by 
taking into account the quartic order terms in the expansion: 


H - -2 cos k—V cos Qx 

^ -2-V + k^ + yeV/2 - ^ {2k‘^ + (Bl) 

where we use the one-dimensional representation and set f = 1 
and the origin of the coordinate at the potential minimum for 
simplicity. 

Let’s define the bosonic ladder operators \n) = 

Vn + 1 |n + 1), a \n) - ^fn |« - 1) and a^a \n) - n \n), with 
which: 


k 

X 


—p[a' - a) 

(fll- + a) 


(B2) 


the canonical commutation relations are preserved. Here tu = 

eV^. 

As a consistency check, we first consider the quadratic or¬ 
der: 


® + (a' -Ha) -(a'-a) 

= — ^a^a + aa^) = cu («-H 1/2) (B3) 


which is diagonal in the \n) representation and has equally 
spaced energy levels. This is the harmonic oscillator. 

We can generalize the expression in terms of the ladder op¬ 
erators to the quartic order: 




(4) _ 




24 

'96 

'96 


1 1 


V 


—I — I a'a ''a^a' — (4« 4 - 6)1 -I a^a ' 4 - H.c. 


1 1 


V 




Jl 1 \ 

4 - 

6 n^ 4 - 6 n 4 - 3 

*(2 ^ y) 


(B4) 


which now has off-diagonal components. Still, for the lowest 
energy levels we can get a sound approximation with a reason¬ 
able up-limit for n and exact diagonalize the Hamiltonian in 
the \n) representation. The results for V = 1.6 and Q — 1/31 
are shown in the blue marks in Fig. |3 

It is straightforward to generalize to higher order expan¬ 
sions. 


Appendix B: Beyond the free-electron approximation near the 
hand edges 

In Sec. m we have shown that by expanding the energy 
dispersion to the quadratic order, the free-electron approxi- 


Appendix C: The tunneling amplitude from WKB 
approximation 

To make a more quantitative estimate of the tunneling am¬ 
plitude between the neighboring wells feff in a slowly varying 
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potential V cos(Qx), we generalize the canonical WKB ap¬ 
proximation to the current dispersion relation: 

2fcos[^(x)] - VcosiQx) - E 
f H-= -ycos(ex)-£ (Cl) 

= -EI2t-(VI2t)cos(Qx) 

+ yl[E/2t + iV/2t)cos(Qx)f-l 

where the energy E is near the lower band edge and x 6 
[0,27r/2] labels the sites between the neighboring wells. In¬ 
side the barrier, the wavefunction decays as where k - ik 
is the imaginary momentum. The tunneling amplitude across 
all sites among the barrier is 

4ff oc exp-^log{-£'/2f (C2) 

X X 

- (y/2f)cos(2x) H- ^[E/2t + (V/2t) cos iQx)f - l| 

Replacing the summation with the integral dx, it is 
straightforward to see that feft exp [-(27rQ:/0] is exponen¬ 
tially suppressed by the distance 27t/Q between the potential 
wells, where 

O' = log I ^{EI2t + (y/2f) cos (y)f - 1 

-EI2t - (y/2f) cos (y)) dy (C3) 

is an increasing function of V/2t, since a higher barrier tends 
to reduce the tunneling amplitude. At the band bottom E = 
-2t - y, the dependence of feff on y/2f and Q is qualitatively 
consistent with Eq. [T5]from simple variational analysis. 


Appendix D: Application to metal-insulator transition of 
generalized Harper’s equation 

It is straightforward to generalize Harper’s equation to in¬ 
clude further neighbor hoppings and higher harmonic compo¬ 
nents in the incommensurate potential: 

t^an'-nfn' + Vfn^b,„ exp {iQm ■ n) = Ef„ (Dl) 

n' m 

where a„ and b„, are given parameters and Q is an incom¬ 
mensurate wave vector. In this appendix, we briefly discuss 
locating the metal-insulator transition separating the insulat¬ 
ing phase for large V and metallic phase for small V in the 
Q/27: <& 1 limit from a two-dimensional Hofstadter problem 
perspective. 

The generalized Harper’s equation in Eq. ID II is equiv¬ 
alent to a two-dimensional lattice problem with dispersion 
f(^) = tYjCin exp {ik^n) + VYjbm exp (ikying in a magnetic field 

n m 

Bz - 2/27r. According to Sec. |II] the states at energy E are 
fully de-localized (localized) if the Eermi surface defined by 
e{k) = £ is open along the y direction (x direction). Yet, as 


we have shown in Sec. [1111 the true metal-insulator transition 
is determined by the quantum tunnelings across the x and y 
directions when the Eermi surface is closed: 



(D2) 


where k^iky) and ky{kx) are solutions of the dispersion e{k) = 
E for a given ky and k^, respectively. It is also straightforward 
to see that ty and are related to the 1/Q coefficients in the 
one-dimensional tunneling amplitude between the neighbor¬ 
ing potential wells in the original Harper’s equation and dual 
Harper’s equation, respectively. We have listed some exam¬ 
ples as follows: 

Generalized Harper’s equations with implicit self-duality: 
In Ref. mobility edges have been identified for the follow¬ 
ing generalized Harper’s equations: 


t{Un-\ + U 

V„=2V 


«+i) "t" — EUfi 

cos (Qn + (p) 

1 - r cos {Qn -H f) 


(D3) 


which can be mapped to an equivalent non-interacting two- 
dimensional lattice model in the presence of an incommensu¬ 
rate magnetic field of - QI2 ji with dispersion: 


2t cos kx + 2V - 

1 - r cos ky 

2t cos ky(l - r cos ky'^ -i-2V cos ky 
—2rt cos ky cos ky + 2t cos k^ 


= E (D4) 

= E{\-rco^k^ 

{2V -H rE) cos ky - E 


It is straightforward to see that for 2V +rE > 2t{2V + rE < 
2t) the constant energy contour at E is elongated along the x 
direction (y direction), 4 (Jy ) is dominant and the model is 
localized (de-localized) along the initial x direction. Clearly, 
there is a mobility edge at £ - 2{t — V) j r, where the Eermi 
surface is symmetric under kx <-> ky and £ = ty. This result 
is fully consistent with the conclusions in Ref. It is also 
simple to check the consistency of our argument with the other 
model in Ref. |9|: 


V„ = 2y [1 - cos {Qn -H f)] / [1 H- r cos {Qn -H f)] (D5) 


Similar self-duality and mobility edge arguments can be 
generalized to two-dimensional incommensurate Hofstadter 
problems with dispersion of the following form: 

a + b ■ h (cosky) 

h (cos kx) H---2-^ = £ (D6) 

c h (cos ky j 

where h{t) is a power series of t and a, b, c are parameters. 
The mobility edge is at £ - b — c. 

Generalized Harper’s equation with nearest and next- 
nearest neighbor hopping: we don’t always have the luxury 
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FIG. 11: Phase diagram and mobility edge of generalized Harper’s 
equation in Eq. |D7] for various next-nearest and nearest neighbor 
hopping amplitude ratio r. The results are obtained through calcula¬ 
tions of Eq. ID2I The solid lines are the exact mobility edges^ for the 
approximants with hopping amplitudes t„ = 


of symmetry for the determination of the relative strength of 
ty and tx in a generic dispersion e(k), therefore, we need to 
calculate ty and tx using Eq. ID2I One simple example is: 

t + M„+i) -H r (m „_2 + Mn+ 2 )] + Vcos(Qn)u„ = £'m„(D7) 

where r is the ratio of the next-nearest neighbor to the near¬ 
est neighbor hopping amplitude. We numerically calculate ty 
and tx from the dispersion relation of the corresponding two- 
dimensional Hofstadter problem e(k) - t (cos kx + r cos 2kx) + 
V cos ky, and the metal-insulator phase boundaries for various 
r is shown in Fig. [TT] In addition, the exact locations of the 
mobility edges^ for approximants of Eq. ID7I with hopping 
amplitudes u = 1,2, • • • are shown for comparative 

purposes. 


Appendix E: Characteristic correlations of a real-space Fermi 
surface 

From the duality transformation in Eq. [3]we obtain: 

{^clcoj - — ^ c^„Cm'&y^p{imnQ + Wm - i6m') (El) 

m,m'=l 


where c is the electron operator for the dual Harper’s equa¬ 
tion with V > 2t. We take advantage of the fact that for a 
quasiperiodic system the correlations are independent of the 
choices of 6: 


{clco) = ^ ^ cl,c„f f dOexpii, 
= y ^ exp (imnQ) 6 „, 


mnQ -t- iOm — Wm') 


= exp {imnQ) = p (nQ) (E2) 

m 

where p{nQ) is the Fourier transform of the electron density 
of the dual model. 


In Sec. HI] we have shown that for Harper’s equation with 
V <2t and Q <& 2n, the electron states near the band edge are 
effectively localized in the physical limit with finite energy 
resolution 6, and the Green’s function ^c|co^ is exponentially 
suppressed as a function of «; otherwise, near the band center 
between the mobility edges, (clco^ is finite for all n. Accord¬ 
ing to Eq. IE2I for the dual Harper’s equation with V > 2t, the 
Fourier transform of the density p (nQ) near the band edge is 
order ~ (9(1) when n is small, since the electron density has 
an induced charge density wave with wave vector (9; how¬ 
ever, at large u —> 00 , p (nQ) —> 0 is suppressed exponentially. 
In contrast, around the band center within the energy range 
E 6 (2f- V, V—2t), the emergence of real-space electron den¬ 
sity singularities (see our conclusions in Sec. lIVl i allows the 
Fourier components p{nQ) to be finite for all n 6 Z. 





